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Abstract 

We report on a search for the decay — > Dfp~ in a sample of 90 x 10^ T{4S) decays into B meson 
pairs collected between 1999 and 2001 with the BABAR detector at the PEP-II asymmetric-energy 
e~^e~ collider. No significant excess of signal events above the expected background is observed. 
We set a 90% C.L. limit on the branching fraction B{B^ Dfp~)< 1.9 x 10^^. Assuming a flavor 
SU(3) symmetry relation between the decays B^ — > Dfp^ and B^ D^p^ , we set a limit on the 
ratio of CKM-suppressed to CKM-favored amplitudes r{Dp) < 9.5 x 10^'^ at 90% C.L. All results 
are preliminary. 
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1 Introduction 



The Cabibbo-Kobayashi-Maskawa (CKM) quark flavor-mixing matrix provides an elegant ex- 
planation of the origin of CP violation within the Standard Model. CP violation manifests itself 
as a non-zero area of the unitarity triangle 2 . While it is sufficient to measure one of the angles 
to demonstrate the existence of CP violation, the unitarity triangle needs to be over-constrained 
by experimental measurements in order to demonstrate that the CKM mechanism is the correct 
explanation of this phenomenon. One of the important measurements is constraining the angle 
7 = avg(—VudV*i^/VcdV*fj) of the unitarity triangle. A measurement of sin(2/? + 7) can be ob- 
tained from the study of the time evolution of the B^^D^*^^^^ and B^^D^*^~ p'^ |1] decays, a 
large sample of which is already available at the B-factories, and of the corresponding CKM sup- 
pressed modes B^^D^*'>^t:~ and B^^D^*^~^ p~ [3]. The first measurements of sin(2/3 + 7) using 
50^^{*)T7r± decays have been recently published 013) and a similar analysis using B^^D^ 
decays is being reported at this conference 0. 

The measurement of sin(2/? + 7) using B^^D^p^ decays requires the knowledge of the ratio 
of the decay amplitudes, r{Dp) = \A{B^^D^ p^)/A{B^^D~ p'^)\. Unfortunately, the direct mea- 
surement of the branching fraction B{B^^D~^ p~) is not possible with the currently available data 
sample due to the presence of the copious background from B^—^D^p~. However, assuming SU(3) 
flavor symmetry, r{Dp) can be related to the branching fraction of the decay B^ Dfp~ [3]: 



r{Dp) 



(tane'c)^! 
IDs 



IB{B^ D+p- 
B{B^ D-p^ 



(1) 



where 6c is the Cabibbo angle, and /d//ds is the ratio of D and Dg meson decay constants [H]. 
Other SU(3)-breaking effects are typically assumed to be of order 30% 6^. 

Since D^p^ has four different quark flavors in the final state, only a single amplitude contributes 
to the decay. The presence of the meson makes such decays easy to identify. Fig. ^ shows the 
dominant Feynman diagrams for the decays B^ D~p~^, B^ — > D~^p~, and B^ — > D^p^. Eq. (|H) 
assumes that the color-suppressed direct ly-exchange amplitude for B^ —>■ D^p~ is negligibly small, 
which is supported by the data 






Figure 1: Dominant Feynman diagrams for CKM- favored decay B^ ^ D p^ (a), doubly CKM- 
suppressed decay B^ — > D^p" (b), and the SU(3) fiavor symmetry related decay B^ Dfp~ 
(c). 



The present limit on the branching fraction B{B^ Dfp-) is 7 x 10"^ at 90% C.L. (TUl. The 
related decay B^ Dfir" has recently been observed [2]. 
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2 Data Sample, Detector, and Simulation 



We use a sample of 90 x 10^ r{4S) decays into BB pairs collected in the years 1999-2002 with 
the BABAR detector at the PEP-II asymmetric-energy S-factory jllj . Since the BABAR detector 
is described in detail elsewhere |12| . only the components of the detector crucial to this analysis 
are summarized below. Charged particle tracking is provided by a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH). For charged-particle identification, ionization energy 
loss (dE/dx) in the DCH and SVT, and Cherenkov radiation detected in a ring-imaging device 
(DIRC) are used. Photons and neutral pions are identified and measured using the electromagnetic 
calorimeter, which comprises 6580 thallium-doped Csl crystals. These systems are mounted inside a 
1.5 T solenoidal superconducting magnet. We use the GEANT [13^ software to simulate interactions 
of particles traversing the BABAR detector, taking into account the varying detector conditions and 
beam backgrounds. 

3 Analysis 

We select events with a minimum of three reconstructed charged tracks and a total measured energy 
greater than 4.5 GeV, as determined using all charged tracks and neutral clusters with energy above 
30 MeV. In order to reject e+e" — > qq, {q = u, d, s, c) continuum background, the ratio of the second 
and zeroth order Fox- Wolfram moments must be less than 0.5. 

The selection criteria are optimized to maximize the ratio of the expected number of signal 
events over the square-root of the expected number of background events, S/ \fB. The optimization 
was performed on large samples of simulated signal and qq and BB background events. 

The mesons are reconstructed in the modes Df^(l)7r~^ , K*^K^ and K^{K^)K'^ ^ with 
- and K*^^K-n+. The K° candidates are reconstructed from two oppositely- 
charged tracks with an invariant mass 493 < M^+^- < 503 MeV/ c'^. All other charged tracks in 
the B meson decay are required to originate from the e~^e~ interaction point. Depending on the 
decay mode of the Df mesons, different selections are used to identify charged kaons. In the 

— > (l)'K~^ {(j)^K~^ K~) decay mode, we identify kaons by applying a pion veto with an efficiency 
of 95% for kaons and a 20% pion misidentification. In the K*^K~^ and Df — > K^K~^ modes, 

a tight kaon selection with an efficiency of 85% and 5% pion misidentification probability is required 
for the candidate from the meson. 

The candidates are reconstructed from two oppositely-charged kaons with an invariant mass 
1009 < Mk+k- < 1031 MeV/ c^. The candidates are constructed from the K' and a 7r+ 
candidates and are required to have an invariant mass in the range 862 < Mj^-^+ < 922 MeV/c^. 

The polarizations of the K*^ and cj) mesons in the decays are also utilized to reject back- 
grounds through the use of the helicity angles 6h(K*^) and The helicity angle is defined 
as the angle between one of the decay products of the K*^ (or (p) and the direction of flight of the 
meson itself, in the meson rest frame. Background events are distributed uniformly in cos Oh since 
they originate from random combinations, while signal events are distributed as cos^ Oh- Both cj) 
and K*^ candidates are required to have | cos Oh\ > 0.55. Finally, after constraining the decay 
products to the same geometric vertex with a probability greater than 10~^, the candidates 
are required to have an invariant mass within 8 or 9 MeV/c^ of the known value |10| . depending 
on the mode. 

The neutral pion candidates are reconstructed from a pair of photons with a minimum energy 
of 30 MeV. The invariant mass of the photon pair is required to be within a window 100 < m^^ < 



9 



160 MeV/c^. After the mass of the vr^ candidate is constrained to 135 MeV/c^, it is combined 
with a track originating from the interaction point to form p~ candidates. The charged tracks are 
required to pass a loose pion selection. We require that the invariant mass of the two pions forming 
the p~ candidate be within 160 MeV/c^ of the known value |10j . 

We also take advantage of the p~ polarization in the — > p~ decays, requiring that the 
cosine of the helicity angle Oh{p) be either larger than 0.3, 0.35, 0.2, or smaller than —0.15, —0.1, 
—0.4 for the Df K*^K^ , and Df — > K^K^ modes, respectively. Signal events are 

distributed as cos^ 6h{p), modulated by the energy dependence of the vr'' efficiency. The asymmetric 
selection takes into account the larger probability to find a random combination of charged and 
neutral pions in the forward direction in 6h{p), which corresponds to low vr" energy. 

We combine oppositely-charged and candidates to form B'^ —>■ Dfp^ candidates. The 
mass of candidates is constrained to the known value [lOj . In order to reject events where the 
comes from a B and the p~ from the other B, we require that these two candidates have a 
vertex fit probability greater than 0.6% of originating from a common vertex. 

We suppress combinatorial background from qq production using the event topology, computing 
the angle {9t) between the thrust axis of the B meson decay product candidates and the thrust axis 
of all the other particles in the event. In the center-of-mass frame (CM), BB pairs are produced 
approximately at rest and produce a uniform cos 9t distribution. In contrast, qq pairs are produced 
back-to-back in the CM frame, which results in a |cos0t| distribution peaking at 1. Depending 
on the background level of each mode, | cos 9t\ is required to be smaller than a value which ranges 
between 0.6 and 0.85. 

We also exploit the characteristic sin^ 6b angular distribution for e^e^ BB events, where 
6b the angle between the B candidate flight direction and the direction of the incident electron 
beam in CM frame. In contrast, the qq backgrounds tend to maintain 1 + cos^ 6b distribution 
characteristic of spin-1/2 particles. We require that | cos of the B candidate is less than 0.6 for 
£>+ ^ (^7r+, 0.8 for L>+ ^ K*^K+, and 0.7 for D+ K^K+. 

We further suppress backgrounds using a Fisher discriminant, J^, constructed from the scalar 
sum of the CM momenta of all tracks and photons (excluding the B candidate decay products) 
flowing into 9 concentric cones centered on the thrust axis of the B candidate ^5]. The more 
spherical the event, the lower the value of Figure 121 shows a plot of the Fisher discriminant for 
signal and continuum background events from Monte Carlo. 

We extract the signal using two kinematic variables rriEs and AE. The first is the beam-energy- 
substituted mass mEs = \J {s/2 + -pi ■ pB)'^/Ef — pg, where ^/s is the total e~^e~ center-of-mass 

energy, (£'j,pj) is the four- momentum of the initial e~^e~ system and ps is the B^ candidate 
momentum, both measured in the laboratory frame. The second variable is AE = E^ — -v/s/2, 
where E'^ is the B^ candidate energy in the CM frame. For signal events, rriES peaks at the B 
meson mass with a resolution of about 3 MeV/c^ and AE peaks near zero with a resolution of 
about 30 MeV, indicating that the candidate system of particles has a total energy consistent with 
the beam energy in the CM frame. For the purposes of selection optimization, we define the signal 
region between 5.27 < tties < 5.29 GeV/c^ and |A^| < 45 MeV. Fi gure El shows the tties and AE 
distributions for simulated signal events and background from decays B^ D*'^p~ reconstructed 
as B^ Dfp-, Df (j)7r+. 

Approximately 40% of the selected events contained two or more B^ Df p^ candidates that 
satisfy the criteria listed above. In such events we select a single B^ candidate based on: 1) the 
reconstructed mass of the Df meson, 2) the reconstructed mass of the vr*^ candidate, and 3) the AE 
variable. The choice is made using these variables in a hierarchical manner with the AE selection 
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Figure 2: Fisher discriminant for simulated signal (solid red circles) events and continuum back- 
ground (open black circles) reconstructed as Dfp^,Df (pii^ 

applied last, in order to avoid creating a bias in the A£' distribution of the background events. No 
significant bias is observed in the large Monte Carlo sample of generic qq production events and 
and B^ decays, excluding h — > ucs transitions. 

After the above selection, two main classes of backgrounds remain. First, there is combinatorial 
background from qq production and generic B meson decays. We describe this background by 
a two-dimensional distribution function of mEs and /S.E. In msS) this background is described 
by an ARGUS function, dN/dx oc x^Jl — 2x^/s exp [— ^ (1 — 2x^ /s)\, characterized by the shape 
parameter ^ ■ In AE', the combinatorial background is well described by a first-order polynomial. 

Second, B meson decays with a Df and a light meson in the final state that correspond to 
b — > ucs quark transitions produce distributions similar to signal events in mES> but are typically 
shifted from zero or broadened in /S.E. These types of background will, hereafter, be referred to 
as peaking backgrounds. We have found that the largest contributions in the signal region come 
from decays B^ Dfp^, B^ D*~^p~ , and B^ D*^tt^ . Of these decay modes, only the latter 
has been previously measured 9 . We determined the mES and AE distributions of these decays, 
reconstructed as signal B^ — > Dfp~ chain, from a large Monte Carlo sample, equivalent to several 
times our data luminosity. 

Our kinematic selection suppresses higher resonances such as /9(1450) and non-resonant B^ — > 
Df'K~7r^ component, which are found to be negligibly small in B^ D^tt^tt^ decays 2i- We 
ignore any possible contributions to DfTr~'ir^ final state other than tt^tt^. 

4 Yield Extraction 

Figure0]shows the distribution of events in the (mES; AE) plane for each of the decay modes. In 
the signal region, we observe 1 event for Df (pTT+, 4 for Df K*^K+, and 2 for D+ K^K+, 
with comparable numbers expected from the backgrounds. 

To extract the signal yield, we perform an unbinned extended maximum-likelihood fit to the 
(rn-ES) Ai?) distributions of all three D'^ decay modes simultaneously in the ranges 5.2 < rn-ES < 
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Figure 3: tues (left) and AE (right) distributions for simulated signal (solid red circles) and ■ 
D*~^p~ (open blue circles) peaking background events reconstructed as Dfp^, Df cpir^ 



5.3 GeV/c^, -200 < < 200 MeV. Tlie total number of events in the sample is 163. The 
likelihood function contains the contributions from the combinatorial backgrounds, signal, and 
peaking backgrounds. The probability density functions for combinatorial and peaking backgrounds 
are found to be common to all three modes in Monte Carlo simulation, while the shapes of 
the signal distributions are determined independently for each Dj" decay mode. No significant 
correlation between niES and AE is observed in Monte Carlo samples, and the likelihoods used in 
the fit ignore any such correlation. 

Eight free parameters constrained by the fit include the shape of the combinatorial background, 
characterized by the parameter ^ in ttt-es and the constant term and a linear slope in AE, as well 
as the combinatorial background yields in each mode and the branching fractions of decays 
B^ Dfp^ and B^ D*'^p~ . The signal and peaking background efficiencies are constrained to 
the values determined by simulation for each decay mode. The branching fraction of B~^ — > 
Dfp^ is constrained to be half of i3(i?'' — > Dfp~) from isospin symmetry, and the branching fraction 
of B^ — > D*'^7r~ , is fixed in the fit to the value measured by BABAR |^. The uncertainties due to 
this assumption are included in the systematics. The results of the fit are shown in Fig. El 

The reconstruction efficiencies, as estimated from the Monte Carlo simulation, and the number 
of observed events in the signal region are summarized in the top portion of Table H The bottom 
portion of Table ^ summarizes the signal, combinatorial background, and peaking background 
yields from the fit to data. The central value of the branching fraction returned by the fit is 
B{B^ D+p-) = [0.2 ± 0.7(stat.)] x 10"^. The log-likelihood of the fit is consistent with Monte 
Carlo expectations. 



5 Systematic Uncertainties 

The systematic uncertainties are dominated by the potential contributions from additional peak- 
ing backgrounds, such as decays S° ^ D*-K+, B+ D*+ p^ , B^ DjK*+, B° D+a{, and 
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Figure 4: /S.E vs mEs distribution for Dfp^ candidates reconstructed in the sample of 90 x 

10^ T{AS) decays in D+ <f>7r+ (full circles), D+ K*^K+ (empty diamonds), and K^K+ 
(black stars) modes. The box corresponds to the signal region 5.27 < mEs < 5.29 GeV/c^ and 
lA^I < 45 MeV. 




Beam-Energy Substituted B Mass (GeV) A E (GeV) 

Figure 5: tues (left) and AE (right) distributions for — > Dfp" candidates reconstructed in 
the sample of 90 x 10^ T(4S) decays. The solid curve corresponds to the full probability density 
function from the combined fit to all Df decay modes. The dashed line is the contribution from 
the combinatorial backgrounds. 
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Table 1: The reconstruction efficiency e, number of candidates A^obs observed in the signal re- 
gion 5.27 < mES < 5.29 GeW /(? and |Ai?| < 45 MeV, the contributions from the signal mode 
— > D'^p~ {Nsig), and combinatorial (A'comb) and peaking backgrounds from the fit to data. 





Dt 0vr+ 


Df K*^K+ 


Dt K^K+ 


All Dt Modes 


e 






3.4% 


1.4% 


2.1% 




A^obs 






1 


4 


2 


7 


iVsig 






0.16 ±0.60 


0.06 ±0.21 


0.11 ±0.41 


0.32 ±0.75 


A^comb 






1.64 ±0.26 


2.88 ±0.32 


0.74 ±0.18 


5.27 ±0.44 




> D*^ 


P 


0.04 ±0.14 


0.02 ±0.07 


0.03 ±0.10 


0.10 ±0.19 




> D*^ 

^ s 




0.42 ±0.29 


0.15 ±0.10 


0.30 ±0.21 


0.87 ±0.37 








0.01 ±0.01 


0.00 ± 0.00 


0.01 ±0.01 


0.02 ± 0.01 



B" — > D*~^a^ . Their contributions are estimated from Monte Carlo. In addition, the branching 
fraction of B^ D*^'ir~, fixed in the fit, is varied within its experimental uncertainties [S]. Al- 
together, the uncertainties in the peaking background contributions amount for a ±0.2 x 10~^ 
systematic uncertainty on B{B^ Dt p^)- Other systematic uncertainties include an 11% relative 
uncertainty in B{Dt (j)7r~^) ^7j) ^ relative uncertainty in reconstruction efficiency due to 
Monte Carlo statistics, and uncertainties in charged track reconstruction efficiency, tt", Kg, and 
charged kaon identification. 

6 Preliminary Results 

Combining all systematic uncertainties, the central value of B^ Dtp^ branching fraction is 

B{B° Dtp') = [0.2 ± 0.7 (stat.) ± 0.2 (syst.)] x 10"^ , 

consistent with zero within the current level of precision. We set a 90% Bayesian confidence limit 
at 

^ Dtp') < 1.9 X 10-^ 

assuming a constant prior for B{B^ Dtp') > 0. The likelihood distribution of B{B^ Dtp') 
and the 90% C.L. limit are shown in Fig. 1^1 

Using Eq. ^ and the recent lattice QCD value fos/ fo = 1-22 ± 0.04 |H1, and assuming no 
additional flavor SU(3) violation, we compute the value and 90% confidence limit on r{Dp) to be 

r{Dp) = 0.003 ±0.006 (stat.) ±0.002 (syst.) 
r{Dp) < 9.5 X 10"^ (at 90% C.L.) 

The low value of r{Dp) compared to r(Z)^*V) determined from B^ — > Di*^^iT' decays jSl IHl 
is somewhat unexpected. It implies small sensitivity of CP asymmetries in B^^D^ p^ decays to 
sin(2/3 ± 7), making that measurement significantly more challenging. 
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Figure 6: Likelihood distribution for B{B^ Dfp^), combining both statistical and systematic 
uncertainties. The dashed lines show the Bayesian 90% confidence level limit assuming a constant 
prior for B{B^ D+p-) > 0. 
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